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Abstract 

The structure of the Roper resonance as a hybrid baryon is investigated through studying the transitional 
ampUtudes in J/ip pN* , N* N* decays. We begin with perturbative QCD to describe the dynamical process 
for the J/tp — * 3q + 3q decay to the lowest order of Os, and by extending the modified quark creation model to 
the J/tp energy region to describe the J/ip 3q + 3q + g process. The non-perturbative effects are incorporated 
by a simple quark model of baryons to evaluate the angular distribution parameters and decay widths for the 
processes J/tp pN* , N* N* . From fitting the decay width of J/ip — > "fpp to the experimental data, we extract the 
quark-pair creation strength gi — 15.40 GeV. Our numerical results for J/tp ^ pN* , N* N* decays show that the 
branching ratios for these decays are quite different if the Roper resonance is assumed to be a common 3q state or a 
pure hybrid state. For testing its mixing properties, we present a scheme to construct the Roper wave function by 
mixing \qqqg) state with a normal \qqq, 2s) state. Under this picture, the ratios of the decay widths to that of the 
J/V' — > pp decay are re-evaluated versus the mixing parameter. A test of the hybrid nature of the Roper resonance 
in J/t/i decays is discussed. 

PACS number(s): 14.20.Gk,13.25.Gv,12.39.Mk,24.85.4-p,12.38.Bx 
Keywords: Roper resonance structure;Hybrid sta.te;J /ip 



1 Introduction 

The theoretical studies on the hadronic structure remain as long outstanding problems in particle physics. About 
four decades ago, the naive quark model was proposed to count for SU(3) symmetries in hadronic spectroscopies[l]. 
After that, much progresses have been made in many theoretical aspects, such as the predictions on the mass 
spectrum, decay rate and electromagnetic coupling properties for hadrons (for details,see recent review[2]). So far, 
much more phenomenological models based on the quark constituents have been developed with a large number 
of successes for further understanding the relations between quarks' degrees of freedom and fundamental field 
of quantum chromodynamics (QCD) theory. Despite these successes the quark models acquired, they are still 
confronted with some challenges. One of the outstanding questions is the fact that, in addition to the conventional 
quark states, QCD theory also predicts the existence of hadrons with 'excited glue', which is called hybrid states 



historically. Therefore, searching for these new hadronic states is of great importance both to QCD theory and the 
common quark models. If no hybrid meson/baryon states are found, it would imply that our current understanding 
of QCD theory should be modified, and the dynamics within quark model would have to be changed. On the other 
hand, if an unconventional baryon state is discovered in experiment, it would be necessary to determine whether it 
is a hybrid state or a common quark state, or their mixture. However, the answer to this question seems to be more 
difficult from experimental aspects duo to the fact that it depends on a comparison with theoretical predictions on 
the specific properties for the new states. 

Therefore, the identification of hybrid states would require a close collaboration between experimental and 
theoretical experts[3]. As our knowledge about common hadronic structures is mainly acquired from their mass 
spectrum, hadronic spectroscopy will continue to be a key tool in this field. Various models and methods have 
been used to predict the spectrum of hybrid mesons/baryons, such as the bag model, QCD sum rule,the flux tube 
model,and so on. Though each model assumes a particular description of excited glue, fortunately they often 
reach similar conclusion regarding the quantum numbers and approximate masses of these states. For instance, the 
predictions on the light hybrid mesons are in good agreement with each other, with the so-called exotic number 
mass about 1.5 ~ 2.0GeV. Experimentally we now have some candidates for J^*^ exotic 
mesons, for example, 7ri(1400) with 1 seen in ?77r[4], and 7ri(1600) seen in pTT,ri-K and 6i7r[5]. In baryon sector, 
the Roper resonance N*(1440) has been suggested to be a potential candidate of hybrid baryons for a long time. 
In the pioneering work of the evaluation of the hybrid baryon mass, Barnes and Close once used the bag model 
to predict that the fightest hybrid baryon was an "extra" ^^N'lPn) state with energy at about 1.6 GeV[6].The 
subsequent identification of this hybrid baryon with the Roper resonance was confirmed basically by Golowich,Haqq 
and Karl[7]. The estimation of its mass was also carried out by using the QCD sum rule, a general prediction came 
with a conclusion that the lightest |^ hybrid baryon is near 2.1 GeV[8]. Recently, this value has been corrected 
by Kisslinger to be about 1.5 GeV, very suggestive of the the Roper resonance [9]. It seems that one might find 
the state with excited gluonic degrees of freedom by studying baryonic mass spectrum. Unfortunately to identify 
a hybrid baryon by the mass spectrum remains an unconquerable difficulty, mainly due to the lack of exotics 
(this situation is quite different from that in hybrid mesons). Hence, for their identification we have to predict 
other experimental observables characterized by their distinct properties as a signature of hybrid baryons. For 
example, strong decay amplitudes[10], photoproduction and electroproduction(EM) amplitudes and so on, can be 
served for this purpose [11]. Especially,the study of the A''* spectroscopy from the J/V' hadronic decays at Beijing 
Electronic-Positron Collider (BEPC) may be an alternative tool to identify the nature of the Roper resonance[12]. 
One might expect hybrid baryons to have a larger production amplitudes from J/'ip decays for hybrid baryons than 
the conventional qqq baryons, because the three gluonic intermediate states produced in J/tp annihilations may 
have a large overlap with final hybrid baxyons. Furthermore, this approach has an additional advantage that it is 
an isospin 1=1/2 filter,so that no A (and hybrid A)states are present to complicate the analysis. 

Recently, BESII has finished datar-taking for 58 million more J/ip events, which is about two order of magnitude 
more statistics than MARKII data, and one order of magnitude more statistics than BESI data. With such statistics, 
partial wave analysis of relevant channels are possible, the angular distributions and the decay widths of processes 
J/tp — »■ piV*(1440) and J/tp — » iV*(1440)A''*(1440) are expected to be available from partial wave analysis in the 
near future. 

What properties are able to serve as a signature to distinguish the Roper resonance between a hybrid baryon 
and a conventional three quark baryon? As well-known, the spatial distribution of constituent components bound 
in the Roper resonance is described as the first radially excited state in the naive quark model, while as a ground 
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may provide us an effective method to test the hybrid nature of the Roper resonance through studying transitional 
amplitude involved in a given J/ip decay process. In this paper we pick out some relevant processes, in which 
the J/il) particles decay into a baryon and an anti-baryon pair including Roper resonance to investigate the Roper 
structure as a hybrid state. 

This paper is organized as follows: in the second section, we present our model of the J/il> decay into a hybrid 
baryon and an anti-baryon, and make a choice of wave functions for hybrid baryons. With this model .in section 
3, we formulate the amplitudes of J/^ decays into pN* and N* N* baryon pairs , and its corresponding decay 
widths and angular distributions for these two decay modes. In the fourth section, the determination of the quark 
pair creation strength from the process J/^ —> 'ypp is formulated. Then, we present our main numerical results on 
angular distributions and decay widths versus the mixing parameter for the two decay modes J/tp —^ pN* , N*N* , 
where the structure of the Roper resonance in two different pictures are assumed. We also discuss our results in 
the last section. Some matrix elements of the transitional amplitude are appended at the end of paper. 

2 Model description 

The exclusive decays of the charmonium have been investigated by many authors within perturbative QCD[13]. 
The main dynamical mechanism is simply assumed that the cc quarks annihilate into a minimum number of gluons 
constrained by the charge and parity conservation. The decay of the J/tp particle into a baryon-antibaryon pair 
is currently assumed to proceed via two steps as illustrated by Fig. 1(a). In the first step, the cc pair annihilates 
into three gluons, followed by the materialization of each gluon into a pair of quark-antiquarks. Since the cc quarks 
annihilate only if their mutual separation is less than about l/rric (w-c is the mass of c quark), which is smaller 
than the non-perturbative charmonium radius, and since the average energy of three virtual gluons is of the order 
1-2 GeV which lies in the QCD perturbative region. Thus this step is usually carried out from perturbative QCD 
to the lowest order of Us- Then in the second step the three quarks on the one hand and the three antiquarks on 
the other hand combine to form a baryon and antibaryon, respectively. The nonperturbative dynamics is included 
by the quark wave functions inside the baryons from the quark model. 

The exclusive decays of the J/tp can also provide a new laboratory to study the Roper properties as a hybrid 
state[12]. Since the strong decay of the J/tp particle is a rich-gluon process, one may indeed expect a possible 
situation that a gluon combines with three quarks or antiquarks to form a physical singlet, a candidate of the 
hybrid baryon, denoted by \qqqg) as illustrated by Fig.l.(b). In this decay mode, a quark-antiquark pair with a 
vacuum quantum number J^^ — 0^^ is assumed to create from the QCD vacuum due to the strong quark-gluon 
coupling. For the sake of simplicity , the coupling of the created quark pair is parameterized into an effective 
strength gi within a modified quark pair creation model. Generally speaking, the strength gi may be dependent on 
the model and the energy level. We feel it is reasonable to determine its value from the radiative decay of the J/tp 
into 7PP, which has an almost identical Feynman diagram as shown in Fig. 1(c) as for the hybrid baryon production 
in Fig. 1(b). 

As in naive quark model, the quarks and the gluon bound in hybrid baryons are treated as constituent compo- 
nents with mass ruq and mg, respectively, and their spatial wave function are chosen to be harmonic-oscillator eigen 
wave functions in the center of mass (CM) system of the baryon. The construction of the hybrid wave function 
is completed by combining the quark and the gluon sector into a color-singlet state, which is indeed included in 
the representation of group SU{S)c^ SU(3)c^ SU(3)c^ SU{8)c- Let I, J denote the quantum numbers of to- 
tal isospin and angular momentum for \qqqg) state with In and J3 , the z-projcction components, then the wave 
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Figure 1: Lowest-order diagrams for {a.)J/ip Iqqq) + |<Z<Z<?), (b)J/'!/' \<1Qq) + \<l<l<ig),ic)J/'ip ^pp 
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where | A'^, /3j3)2/,2j is the color-octet qqq part in the hybrid state, with quantum number of the total isospin / 
,z-projection component I3 and angular momentum j ,z-projection J3. a is the color-index, and \g, IJ3 — js)" is 
the color spin wave function for the gluon. 

Especially, for 1=1/2, the wave functions for the color-octet qqq part[6] are 



^8:iiV,/3j3)?,3 = ■^{cb'rc~r^^)x 
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Here (j), x and i/i denote the flavor, spin and color wave functions, respectively, s, A and p denote the symmetry or 
mixed symmetry representations of the total permutation of particles (A is symmetric and p antisymmetric under 
(23) interchange). The study of the electromagnetic transitional properties suggested that the wave function for 
the \qqqg) state should be chosen as a mixing of '*8 and ^8 with an equal ratio[ll],i.e. 



\qqqg) = ^ (| A^g, /3 J3)i,3 -I- \Ng, I^Js)!^!) x symmetric spatial part. 
V2 



(4) 



In the constituent quark model, the wave function of the nucleon and the Roper resonance are assigned as 
a ground state and the first radial excitation of three constituent quarks in the harmonic oscillator potential, 
respectively. The nucleon wave function in the color-spin-flavor space reads. 



|p) = \qqq, Is) = — r=-0''(x''0'' + X^0^) ^ symmetric spatial part, 



(5) 



are the spin and flavor wave functions for the 



^/2 

where is the asymmetric color wave function ,x''(x^) ^.nd 0'' 
qqq cluster with p (or A)-type mixing symmetry, respectively. 

The identiflcation of the Roper resonance as a \qqq,2s) state meets some difficulties in the reproduction of the 
mass spectra and the photoproduction amplitudes. On the other hand, some theoretical studies have predicted 
the existence of a hybrid state at this mass region [6, 8, 9]. So the mixing of those two conflgurations should 
be considered. Li[ll] once mixed the \qqqg) part with \qqq, Is) in the study of A'^*(1440) in the photoproduction 
process, where the wave function of nucleon |p) and the Roper resonance lA'^*) are assumed to be linear combination 



nontrivial order is expressed by \p >= [\qqq, Is) — ^\qqqg}]/-\/l + A low lying hybrid baryon is assumed to be 
the Roper resonance, whose wave function is chosen as to be orthogonal to the nucleon wave function,i.e.|A''* >= 
1*) + \qqq9}] I + In this scheme ,the wave functions for the nucleon and Roper with gluonic freedoms 
are constructed by introducing a phenomenological mixing parameter ^ to combine the gluonic freedom with the 

conventional proton. Its lowest limit ^ = corresponds the extreme situation,namely, the nucleon is the conventional 
|3(/, Is) state, while the Roper is a pure hybrid state. However,the upper limits of the mixing parameter ^ is 
unrestricted. If the fraction of \qqqg) part grows larger in value, it will change the nucleon properties. On the 
other hand, there are two states \qqq,2s) and \qqqg) at the Roper mass region theoretically. These two states may 
mix together due to the residual interaction. One may think that the mixing of the two states might change their 
eigen energies. So we propose an alternative scheme to construct the wave function for the Roper resonance. It is 
assigned as a mixture of the \qqq, 2s) part with the \qqqg) part, 

\N*) = 6\qqq, 2s) + ^/ 1 - 5^\qqqg) , (6) 

which is automatically orthogonal to the conventional nucleon state. Here 5 is a mixing parameter running from 
up to 1. \i 5 = 0, it indicates that AT* (1440) is a pure hybrid state, and 5 = l,the pure \qqq, 2s) state. 

3 Decay widths and angular distributions 

The full decay width of the J/zf; particle into three gluons has been investigated many years ago. Since the 
charmonium decay process involves two different physical scale; the binding energy e of the bound state and the 
heavy c-quark mass mc{mc ^ e), now it is widely accepted that the transitional amplitudes for charmonium decays 
can be expanded to order 0{p/m) , where p is the relative momentum of the cc quarks. For J/ip particles described 
as a Is bound state, only the first order term contributes dominantly to the decay width. This approximation 
results in that the decay width is proportional to the value of the J/ip wave function at the origin |$(0)|. To 
study the exclusive decays of J/ip into a baryon anti-baryon pair, we focus on the dynamical behavior of the gluon 
properties, the outgoing quarks and the information on the baryon structure. We feel that the non-perturbative 
properties of the bound cc quarks and their dynamical effects do not play important roles in the J/tp exclusive 
decays at least in the limit of the non-relativistic approximation, since in which the contribution from cc quarks 
are independent on their relative momenta. Hence the non-perturbative effects within this energy level, together 
with the contribution from the cc quarks and the decay constant of J/ip and so on are all parameterized into an 
overall constant Co. 

The gluons from the cc annihilation play a different role in our model. In the process J/ip —^3g~>3q + 3q, 
each virtual gluon behaves with equal status due to the fact that each gluon finally creates a qq quark pair. While 
in the process J/ip —^3g—^3q + 3q + g, of all three gluons, there being one gluon does not create a qq quark pair 
,and instead is bound inside the hybrid baryon in the final state. We consider this gluon to be a constituent one in 
the bound state. 

If we assume the Roper resonance is a pure hybrid state, i.e. without mixing with \qqq) component in its wave 
functions, there only being one precess as shown in fig. 1(b) should be considered. In general, the Roper resonance 
is assumed to be a mixture of the \qqq} and \qqqg) components. In this case of mixing picture, to study the decay 
of J/ip into the Roper, it is essential to calculate the basic transitional amplitudes of both decay modes given in 
fig.l.(a) and (b). The basic amplitudes for these two processes can be expressed as, 

< iVAr|TiCi|J/V^''^ >.,,4 = < *jv(g',s;)*jv(g,s)|CiTi|J/V'^'') > 
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X <gi,Si,g^,sJ,i = 1..3|CiTi|J/V'(''^ >, (7) 



for fig.l.(a) ,and 



< iViVg|T2C2|J/V'^^ >.„4 = <*h(q',4)*^(g,s.)|C2r2|J/V^''^ > 

j = 1..3 



= X! / *^ n <*/»(?', si)*jv(?,«2:)|gi, Si, = 1-4, j = 1. .3 > 

Si,s^ i=1..4 



X <qi,Si,q'j,s'j,i = l.A,j = 1..3\C2T2\J/ip^''^ >, (8) 

for flg.l.(b). Here N and iVg stand for the qqq and g^gg states, respectively. '^N{q',s'^) and are the 

flavor-spin-spatial wave functions for qqq and qqqg cluster, respectively. 

The hard scattering part in above equations can be written explicitly by the standard Feynman rules, i.e. 

(gi + q'jy(q2 + q'2r{q3 + q^V 
X '«(gi,s'i)7^v(gi,Si)M(g2,S2)7''w(g2,S2)M(g3,S3)7'''y(g3,S3), (9) 

and 



(\i)p (A)A gti.\9up + QvXQtip + gpXQuv 
[qi + q[nq2 + q!2) 
X -"(gi, si)7''«(gi, Sl)M(g2, S2)7''w(g2, S2)M(g3, S3)'f (gs, S3), (10) 



< gi,ai,g^.,s;|C2r2|J/V'(^) > = Co < C2 > {igfgje^'^Kj^^ 



where e'y'^J^, e'-'^*''' arc the polarization vectors for J/ib and the gluon with helicity values A and \i, respectively. 
We perform the calculation in the J/i/; rest system. For J/i/; produced in e^e~ annihilation, its helicity is limited 
to be A = ±1. qi{q'i) and Si(s^) arc the four vector momentum and spin z projection of the anti-quarks (quarks), 
respectively. (Ci) and (C2) are color factors corresponding to the decay mode as depicted in fig.l.(a) and (b), 
and g — ^^A-kUs is a strong coupling constant, gi is an effective strength of the created quark pairs in the ^Pq 
quark model which will be determined from the decay mode J/tl) ^ fpp. Other constants and contributions 
from charmonium bound properties, which are assumed to be independent on a given process, are all put into a 
single overall constant Co- u{q'i,s'i) and v{qi,Si) are assumed to be free Dirac spinors for quarks and anti-quarks, 
respectively. They are explicitly given by 

, , lmi + E'. f 1 \ , Irrii + Ei ( mi+E] \ 

where Ei = qf is the energy of a quark. 

The evaluation of the color factor < Ci > and < C2 > for these decay modes is straightforward in terms of the 
SU(3) structure constant dabc and the SU(3) color matrices T, they are 



^ /-I ^ ^ 1 ijkrjiarpb rrtc Imn. ^ /'10\ 

24 a/3 ^il-^jmJ-kni - vJ-^j 



1 I ijkrjiarpb rpc Imn ^ 

< C2 >= ^dabce'^'TLnmSkirailmn) = ^, (13) 

where T^"'\a = 1..8) are Gellmann matrices, and '4)''{lmn) is the color-octet wave function for three quarks in the 
Igggs) state with antisymmetry by exchanging quark (2) and (3), while the A-type wave function ijj^{lmri), which 



symmetric properties of Gellmann matrices. It is worthy to note that the ratio {C2) / {C\) w 0.2 from Eqs. 
and (13), which means that the decay of the J/ip particle into the hybrid state is a color suppressed process. 
The components involving Dirac spinors in Eqs. (9, 10) can be expressed more explicitly as 
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(14) 



AE^Ei 

where Ei and rrii arc the energy and the mass of the quark. 

From these amphtudes,one may obtain elements of the transitional amplitudes M'^J^]^, for processes J/%1> ^ 
pp,pN* , N* N* by projecting them on the wave functions of the nucleon or the Roper resonance ( see Appendix 
B). In our calculation, we ignore the amplitude for J/ip — » \qqqg) + \qqqg} because we argue that the contribution 
of matrix elements for this process is trivial due to the color factor and an additionally suppressed factor e~^<=/^'' 
in the spatial wave function for the \qqqg) cluster. Then the differential decay width for J/tp — » BB{B : baryon) 
can be expressed as, 



dTjJ/i,'^^'^ BE) 



1 
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(15) 



dQ 327r2 

where Mj/^ is the mass for J/ip, and the if is the momentum of outgoing baryons. 

From the symmetry consideration, the differential decay width of J /xp — > BB can be expressed in a more general 
form, 



drjj/^'-^^ BB) 



dQ 



ascos^ 9), 



(16) 

BB, and can be 



where Ngg is a constant directly related to the experimental total branching ratio of J/t/j 
used to fix the overall constant Co, and 9, is the angle between the positron beam direction and the direction of 
the outgoing baryon. as characterizes the angular distribution of J/'0 — » BB, which can be extracted from the 
experiments. 



4 Strength of quark-pair creation in J/il) decays 

The quark-pair creation model was first introduced by Micu[15] in 1969 in a study of meson decays, which suggested 
that the strong decay of mesons proceeds through a simple quark-antiquark [qq) pair creation from the vacuum, 
with a quantum number J'"' = 0"*"^. With a qq pair creation operator and a dimensionless constants 7, the widths 
of the two-body hadron decays can be evaluated easily from definite transitional amplitudes. Many authors have 
developed Micu's original suggestion and applied it extensively to a number of baryon and meson decays with 
considerable successes. In this model, the constant 7 is a parameter and is obtained from the fitting to the light 
meson or baryon decay widths. 

For the OZl rule allowed decays of the J /xp particle, the qq pair might be created from either a gluon annihilation 
or the QCD vacuum . However, in the energy region of the J/V* decay the strength of the qq pair creation from the 
QCD vacuum may be different from that in the light meson or baryon decays. Therefore, we adopt the modified 
quark-pair creation model to extract the strength of the quark-pair creation by evaluating the decay width of the 



The radiative decay J/xj) — » 7pp is assumed to proceeds via two steps as illustrated by fig. 1(c). In the first step, 
the cc pair annihilates with the emission of a photon and two gluons. In the second step, each one of the gluons 
creates a quark-antiquark pair and another quark-antiquark pair is created from the vacuum, which can be described 
by quark creation model with the quark pair creation strength gi. By a complicated final state interaction, a proton 
and an antiproton are formed. This complicated hadronization process is accomplished by the bound state quark 
wave functions of the nuclcon, which is constructed in the constituent quark model. 

In our quark pair creation model, the quark-antiquark pair with any color and flavor can be created anywhere 
from the QCD vacuum with equal strength. But only those pairs whose color-flavor wave function and space wave 
function overlap with that of baryons in the final state can make a contribution to the decay width. Following the 
usual procedure, the Hamiltonian for the quark pair creation can be defined in the modified ^Po model[16] in terms 
of quark and antiquark creation operators 6"*" and d"*" , 

Hi= ^ j d'kd'k' gi[u{t\s')v{l^,s)]hl0's)d^^0,s)5^{k (17) 

where a(/9) and i(j) are the flavor and color indexes of the created quarks (anti-quarks) , and u{k',s') and 
v{k,s) are free Dirac spinors for quarks and antiquarks, respectively. They are normalized as u^{k,s)u{k,s') = 
v~^(k, s)v{k, s') = 5ss' ■ Ci = 6ij is the color operator for qq and gi is the strength of the decay interaction, which 
will be extracted from fitting the width T{J/tp — » 7pp) to the experimental data. In the non-relativistic limit, gi 
can be related to 7, the strength of the conventional ^Po model, by gi = 2mq7[16]. 
We first define an operator for J/tp f + 2q + 2q transitions, i.e. 

OiJ/i. - J,qq,qq) = CoC{^e)i^g) e^,/^, . {<n + gi^iq. + q',)^ 

■uil'iy s'i)-y''v{qi, si)w(g2, 52)7'' ^^(92, S2), (18) 
where qi{i = 1,2,3) and q'i{i = 1,2,3) are the four-rnornentum vectors for anti-quarks and quarks, respectively. 
Co is an overall constant which can be determined from fitting the decay width of the process J/ip ^ pp. C — 
j^SabTl^n'^jm IS the color operator in terms of SU(3) color matrices T, and p is a strong coupling constant, e'-^*-' 
is a photon polarization vector with helicity Ai — ±1. e j'^^ is the J/tp polarization vector. 

Inside a baryon there are strong interactions among three constituent quarks. To form a baryon from the three 
quarks is a non-perturbative process. Following the usual approach, we account for all strong interaction effects by 
the bound-state quark wave-functions in the final ( anti-)baryon state. The amplitude for J/tp — » jpp transition 
can be written as, 

n'-^\sz,s'^,pi,P2,P3) = {'^p{sz,Pi,qi,q2,q3)tpp{s^,p2,q[,q2,q'3)\OHi\%pf^^{P)), (19) 

where Sz and s'^ are the spin z-projections for the anti-baryon and baryon, respectively. tpp{sz,pi,qi,q2,q3) and 
V'p('Sz)P2, 91,112, 13) are the spin-flavor-spatial wave functions of the baryon and anti-baryon, respectively. In general 
, the structure of spin-flavor wave functions of baryons can be constructed in the constituent quark model. The 
spin and flavor wave functions of the proton and antiproton are similar, i.e. 

where x'' and are the spin-i wave functions of the quark pair with mixed-symmetry (see Eqs. (31,32) in Appendix 
A). 

The spatial wave functions for the proton or antiproton are chosen as simple harmonic-oscillator eigenfunctions 
in their center-of-mass(c.m.) system. However, our calculation is carried out in the J/ip rest system, where the 



function from the cm. system of the baryon to the laboratory system. Here we only perform the Lorentz boosts 
for the spatial wave function and ignore the Melosh rotations of the quark spinors. This approximation will be 
discussed in the next section. 

The decay widths for the process J/V* — » 7PP can be expressed as: 

r = wE E /iMA,....,P.P.P3)rn(^^^(p-p.-p.-P3), (21) 

A Sz,s^/ ,Xi i=l 

where P and M are the four momentum vector and the mass of the J/ip particle, respectively. In the J/ip cm. 
system, P = 0. Ei,E2 and E3 are the energies of the anti-baryon, baryon and photon, respectively. 

The photon in the process J/ip — » 7pp may also come from the bremsstrahlung of the outgoing proton or 
anti-proton following the strong decay J/ip pp. However,this bremsstrahlung mainly contributes to the low- 
energy photons[17]. On the other hand, the existing measurements for the radiative decays are performed with the 
condition rn{pp) < 2.79GeV/c[18], which implies that photons with energies smaller than 130 McV are cut off. In 
our calculation wc ignore the contribution from the bremsstrahlung in the final state to the radiative decay width 
and perform the calculation with the experimental condition to evaluate the quark-pair creation strength. We use 
the branching fraction Br{J/ip 'ypp) = 3.8 x 10~*[19],the light quark mass ruq = 0.22GeV and the harmonic 
constant O.OSGeV^ to extract the strength gi = 15.40GeV. We also change the harmonic oscillator parameter from 
O.OSGeV^ to 0.16GeV^ and find that the strength gi is not sensitive to the choice of this parameter. 

This value is much larger than that from fitting to the light meson decay widths by Ackleh eto/.[16]. Prom 
their calculation, we get gi = 7.3GeV. The difference may indicate that the strength gi is energy-dependent. 
Besides, in Ackleh's calculation, the non-relativistic limit Eq = ruq is used and the Lorentz contraction for the wave 
functions of outgoing mesons is ignored. While in our evaluation, the relativistic transformation on the baryonic 
radial wavefunctions between two different systems are carried out. 

5 Numerical results and discussions 

Although the decay of the J/%p into three gluons and the creation of a quark-antiquark pair from one gluon may be 
treated by perturbative QCD, the formation of a baryon or hybrid baryon is certainly beyond the asymptotic regime. 
We account for this effect by explicitly including the bound state wave function for the three-quark cluster or quark- 
gluon cluster. As in the naive quark model, we treat the quarks and gluon in the baryon as constituent components 
with masses niq and rUg, respectively. We assume that the momentum distribution of the bound components in the 
cluster \qqq} and \qqqg} can be described by wave functions 4>ni{Pp,px) and (l>h{Pa,Pb,Pc) ,respectively. Since the 
spin-independent potential is used in the constituent quark model, the spatial wave functions of the constituents 
can be chosen as simple harmonic oscillator eigenfunctions in their center-of-mass (CM) system [see Eqs. (42,43)]. 

In the CM system of the J/ip, the two baryons are moving in opposite directions with highly relativistic speeds. 
One has to transfer the quark wave function from the CM system of the baryon to the laboratory system. The two 
baryons' radial wave functions in the laboratory system are related to the wave functions in their CM system by 
the Lorentz transformation [20, 21]. In principle the quark spin wave functions have to be changed by the Melosh 
rotations of the quark spinors simultaneously. Por simplicity, this relativistic transformation can be simplified by 
neglecting Melosh rotation effects in the framework of the naive quark model as in the studies of the resonance 
electroproduction[22]. Since the non-relativistic descriptions for constituent quarks are employed in our calculation, 
the momentum dependence of the spin wave function can be ignored. On the other hand, for the J/ip decay process 
we can take the z axis in the direction of the moving baryons. So the third component of the quark spin remains the 



we found that this treatment can be compensated by taking the parameters describing the baryonic properties as 
effective values from fitting to experimental data. Therefor, in our present calculations, we only perform the Lorentz 
boosts for the spatial wave function as follows. 



d{ Up, k' x) 



<t>{kp,kx), (22) 

1/2 

In most quark model calculations, the value of a has been chosen in the range of 0.06 ~ 0.22GeV^, which 
corresponds to the nucleon radius in the range of 0.52 ~ 0.98fm. The range of a value in our model has been 
determined elsewhere by fitting the angular distribution of J/i/> pp[23]. We take the mass of the light quark 
ruq = 0.22GeV and a = O.OSGeV^, which corresponds to the central value of the angular distribution parameter 
Oipp = 0.62 for ,J/%1) pp. The assignment of a mass to the constituent gluon is model-dependent. For instance, 
the lattice QCD calculations predict the mass of the constituent gluon rUg = 0.5 ~ 0.6 GoV . But from the 
phonomcnological theory, the mass of gluons are commonly taken in the range rrig = 0.3 ~ 0.7GeV at the level of 
J/tp decays[24, 26]. In our model, if we assume that the Roper is composed of one constituent gluon and three light 
quarks in Is state, we approximately take the constituent gluon mass as rrig ~ Mn* — Mp w 0.5GeV. In our present 
calculation, the gluon is regarded as a relativistic constituent particle, its mass can also be adjusted to a lower value 
rUg = O.STGeV as taken in Ref.[25]. The coupling constant is chosen to be as = 0.28[14], and the strength of the 
quark pair creation, gi, is determined from the calculation of the decay width T{J/tp — > 7pp) in the previous part. 

5.1 Roper as a pure hybrid bairyon or conventional baryon 

At first, we consider the Roper resonance as a pure hybrid state,namely a qqqg bound state. With selected 
parameters varying within our parameter window, we get the angular distribution parameters a* = 0.42 ~ 0.57 
and a** = (—0.1) — (—0.9) for the decay modes J/tp — » pN*, iV*iV*, respectively. 

Although the dynamical behavior of gluons and created quarks are important to the decay mode, we find that 
the structure of the bound cluster play a dominant role in the evaluation of contributions of amplitudes to those 
decay processes. Varying the mass of a constituent gluon from 0.37GcV up to 0.50GcV, we find that the amplitude 
for the hybrid decay modes is much smaller than that for the normal decay mode, and one obtains 

I < qqq, Is; qqqg, ls\T2C2\J/4>'-''^ > U,^ 



\<qqq,ls;qqq,ls\T-,Ci\J/i!M^ > 
\ < qqq,2s;qqqg,ls\T2C2\J/t''^^ > |, 



< 0.15; (24) 

< 0.04. 



I < qqq, \s;qqq, \s\TiCi\J /ipW > |,^, 

This implies that if the Roper is a pure qqqg hybrid, the decay widths for the process J/ip — > pN* , N*N* are almost 
less than 2% and 0.2% of that for the process J/tp pp, respectively . Where doses this suppression come from? 
A great part is from the depressed color factor due to the different dynamical behaviors of the quarks and gluons 
in these two decay processes. One could imagine the fact that if a gluon to combine with a color-octet three-quark 
cluster to form a color singlet object, not all colored gluons are allowed. This situation is well understood from the 
calculation of the color factor < C2 >, which is almost the fifth of the color factor < Ci >. This indicates that J/ip 
decays into a hybrid baryon is a color depressed process. While in the normal decay of the J/tp particle, the gluons 
and the constituent quarks appear with an equal status due to the fact that each of the three gluons decays into 
a quark and an anti-quark pair and then three quarks (anti-quarks)form a color singlet baryon (anti-baryon). The 
color factor for this process is much larger. In other words, the decay mode for a hybrid baryonic production is color 



depressed process. Other source is coming from the different behaviors of the constituent quarks and gluons bound 
in the baryon. We simply assume that the proton is a pure three-quark system and the N* is a pure qqqg system. 
To calculate the decay amplitude one has to project the basic amplitude onto the bound state wave functions of 
three quarks (and the gluon). The overlap of the basic amplitude with the wave functions of three quarks is much 

larger than that with the wave functions of three quarks and a gluon. By adjusting the gluon mass from mg = 0.37 
GoV up to 0.5GcV, and using the harmonic oscillator parameter a = O.OSGeV^, we finally obtain the ratios of the 
J — > piV* and J/V' — > N*N* decay widths to the width of the process J/ip —^ pp, 



r(j/v 



(A) 



r(j/v,(A) ^ pp) 



r( j/v^(A) ^ pp) 



1.4 X 10"^ {nig = 0.50GcV); 

1.8 X 10''^ {nig = 0.37GeV); 

1.2 X 10"^ {mg = 0.50GeV); 

2.3 X lO""* {rUg = 0.37GeV). 



(25) 



We also change the harmonic oscillator parameter of nucleons a from 0.06GeV^ to 0.21GeV^. The results are 
plotted in Fig. 2. 

Another extreme situation is that the Roper resonance is a pure common 3q state, namely 5 = 1 in Eq.(6), 
which is assigned as the first radial excitation state \qqq,2s}. In this case, we deal only with one decay mode, 
J/^p \qqq) + \qqq) as illustrated in fig. 1(a). With the same parameters for the quark masses and the harmonic 
oscillator the numerical results are 

r(j/^(A) ^pjv*) 



r( j/.0(A) pp) 

r(j/V'<*^ -^N'-N") 



t{j/'4jW pp) 



2.0 ~ 4.5, 
3.2 ~ 22.0, 



(26) 



and their angular distribution parameters are a* = 0.22 ~ 0.70 and a** = 0.06 ~ 0.08, respectively. This indicates 
that the the J/V* — > pN* decay width is much larger if the N* is a normal three-quark excitation than that if the 
A^* is a hybrid. 



5.2 Roper in the hybrid picture mixing with the \qqq,2s) state 

If the Roper resonance is taken as a mixture of a normal \qqq, 2s) state with a pure hybrid state, the decay widths 
for J/ip pN*,N*N* decay modes can be evaluated numerically from Eqs.(39)~(41). The predicted ratios of 
decay widths in logarithmic scale arc plotted in Fig. 2 iu terms of the mixing parameter 5. The dashed, solid and 
dotted curves correspond to the different choice of the harmonic oscillator parameters a = 0.06GeV^, O.OSGeV^ 
and 0.21GeV^, respectively. The light quark masses are taken as niq = 0.22GeV. For each type of curves in small 
5 region, the upper and lower curves correspond the different choice of the gluonic mass rUg = 0.37 or 0.50GeV as 
illustrated in the figure caption. It is clear to see from these results that the ratios grow rapidly with the increase 
of the mixing parameter in small 5 region. While in larger 5 region, the ratios almost remain the same as the 
variation of the gluonic mass from 0.37GeV up to 0.50GeV. This behavior is easily understood from calculations 
of the contributions of the hybrid decay mode to the transition amplitudes. Since the J/V* decays into a hybrid 
baryon are color-depressed processes, the normal decay modes contribute to the transition amplitudes dominantly. 
Hence comparing the calculation of the ratios in J/il) decays with the measurements might provide a new approach 
to probe the structure of the Roper resonance. 







0.5 



1 



1.5 



2 







0.5 



1 



1.5 



2 



5 



8 



(q) 



(b) 



Figure 2: The ratios of decay widths are plotted as a function of the mixing parameter S. The 
A''* (1440) is assumed as a mixture of the \qqq,2s) and \qqqg) states with mixing parameter 
5. The dashed, soUd and dotted curves correspond to the choices of the harmonic oscillator 
parameter a = 0.06GeV^, O.OSGeV^ and 0.21GeV^, respectively. The values for light quark 
masses are taken as rUq = 0.22GeV, (a) The ratio T{J/tjj — > pN*)/T{J/il) — > pp).For each type 
of curves, the upper one corresponds to the mass of the constituent gluons rUg = 0.50GeV and 
the lower one is related to the choice rUg = 0.37GeV . (b) The ratio r(J/V' N*N*)/T{J/ip 
pp). Fov each type of curves, the upper solid curve, dash curve and lower dotted curve correspond 
to the mass of the constituent gluons nig = 0.37GeV and the other curves are related to the 
choice rUg = O.SOGeV . 



6 Conclusion 



Wc have studied the decay widths and angular distributions of the processes J/ip ^ pN* , N* N* in constituent 
quark modeL The strong decay of the J/tp is assume by the emission of three gluons. The transitional amplitude of 
dynamical process for J/ip —> \qqq) + \qqq) is numerically evaluated from perturbative QCD approach by assuming 
that each gluon creates a quark-antiquark pair. While in the process J/ip — > \qqq) + \qqqg),we suppose that two 
pairs of quark-antiquarks are created from the annihilation of two gluons and a light quark pair, uu or dd is created 
from the QCD vacuum with strength gj. Then the three quarks ( anti-quarks ) or three quarks and a gluon are 
respectively combine together to form an ordinary baryon or a hybrid state by some final-state interactions. The 
determination of the strength gi is also carried out with the modified quark creation model by fitting the decay 
width for the process J/ip — » 'ypp to the experimental value, and the numerical result gives the strength gi = 15.40 
GeV by using the same experimental cut-ofT condition on the restriction of the outgoing photon energy. Then the 
transitional amplitudes and angular distribution parameters for the processes J/ip pN* , N* N* are calculated 
numerically. If the Roper resonance is assigned as a pure hybrid state, our numerical results show that the ratio 
r{.J/i>^^^ ^ piV*)/r(J/V'''^' ^ pp) < 2%, and r{J/i;^^'> N*iY*)/r{J/tp^^'> pp) < 0.2%, and their angular 
distribution parameters are a* = 0.42 ~ 0.57 and a** = (—0.1) — (—0.9), respectively. However, when the Roper 
resonance is assumed to be a common 2S state, the results are quite different, with r(J/V'^^^ pN*)/T{J — »• 
pp) = 2.0 ~ 4.5, and r(J/V'^'^^ -» N*N*)/r{J/tp'^'^^ -» pp) = 3.2 ~ 22.0, and with the angular distribution 
parameter a» = 0.22 ~ 0.70, a«« = 0.06 ~ 0.08. This implies that, not only the dynamics of three gluons and 
created quarks , but also the structure of the final cluster state , i.e. \qqq} or \qqqg}, play important roles in the 
evaluation of the amplitudes in these decay processes. So it is suggestive that an accurate measurement of the 
decay widths and angular distributions of these channels may provide us a novel tool to probe the structure of the 
Roper resonance. In our present calculation we also consider the mixing of the normal Roper state with the hybrid 
state. Under this picture, we calculate the decay widths for these decays versus the mixing parameter 5. 

In conclusion, we find that the contributions from the Roper structure to the transitional amplitudes paly an 
important role in evaluating the decay widths and angular distribution parameters for processes J/zp pN* or 
N* N* . Hence, the study of the transitional amplitudes for these decay modes may provide us a new laboratory to 
justify the nature of the Roper resonance. If an accurate experimental data on decay widths and angular distribution 
parameters for these processes are available in the future, one may discriminate the hybrid Roper baryons from 
conventional baryons. 
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10225525, 10055003 and by the Chinese Academy of Sciences under project No. KJCX2-SW-N02. We thank 
E.Gset, J. Hiiefner,H. J. Pirner, Z.Y.Zhang, RN.Shen and J.G.Bian for useful discussions and comments. 



APPENDIX 



A The construction of the hybrid state wave function 

Let I, J denotes the quantum number of total isospin and momentum for \qqqg) state with /a and J3, the z-projection 
component, and the wave function of hybrid state denoted by \Ng,I;iJ[i,)2i2.j can be constructed according Eq.(l) 
by combining the gluonic part and the three quarks part's flavor-spin-color wave function into a color singlet 
state. lierc, we use the dummy index sum as Xya=i 8' Especially, for I=J=l/2,we give explicitly their form in 
spin-flavor-color space as follows; 

6'"' V^^"^ va/A.-1/i:, (27) 

1^^'^"^^" = -i=|l)"-i=(0VS-0>^)X-3/2 
1 1 



-|0)"-i=(<^VS-rC)x-i/2 

l)"-^(0VS-0''V'a)Xv2, (28) 
|iVfl,ii)ii = yf|l)"i[(<^V_v2-<^'x-i/2)^S-(0''X-i/2 + <^V_i/2)v^^] 

- yf|O)"i[(0''X?/2-0V/2)V'S-(0''X^/2+0'x?/2)V'a], (29) 
= y^|O)"^[(0V_l/2-</''X-l/2)V'S-(0^X-l/2+</'V_l/2)V'a] 

I - ^r\Wx'[i2 - 4>Wi2W<. - {4>''xi/2 + 4>\',/M (30) 

where 0, x, V' are respectively responsible for the flavor ,spin,and color wave function of the three quark parts in 
\qqqg) cluster. The index s, A, p denote total permutation symmetry or mixed representation. The mixed symmetric 
states are define as: 

ptype: \xyz)'' = ^{\xyz) - \xzy) + \yxz) - \yzx)), (31) 

A type : \xyz)^ = ^^j={\xyz) -ir\xzy) + \yxz) -ir\yzx) -2\zxy) -2\zyx)), (32) 

if x=y, we must replace \xyz) + \yxz) with \/2\xxz) in the above equation. For example,the mixed symmetric spin 
wave functions read: 

x?/2.i/2 = ^(im>-inT)), (33) 

Xv2,i/2 = ^(1 TTi) + I nt) - 2| UT», etc. (34) 

B The matrix elements of the amphtude for J/'^ \qqq) + \qQQ9) 



Since the hybrid wave function is chosen as a mixture of the *8 state with the ^8 state with equal ratio,e.g. 



J/V* — »■ \qgq) + \qqgg)- Note that < tp''\C2\4>^ >= O; where V'" is the asymmetric color wave function of baryon , so 
we only calculate the part < V*" 1*^2 IV''' >• 

+ -^|1 >" {r\Cm{x^\T\X-^} (35) 



< N\C2n\Ng, i - i >i3 = ^|1 >" {r\Cm{x'\T\x'.i) 



+ il - 1 >" {r\Cm){x^\T\x\) (36) 



\/24 



< iv|C2T|ivg, 11 >n = ^ii >" [{ricmix'iTix'..) - {r\C\r.){x'\T\xti}] 

--L^\o)-[{r\cm{x'\T\x'i) - iricmix'iTixi)] 

2v6 2 2 



(37) 



< NianiNg, 1 - 1 >ii = -^1 - 1 >" KV'''iqv'g)(xiTix|> - iricmix^mxp] 

+J^\or[{r\c\K}{x''\T\x\} - {r\cm{x'\T\xt^}] m 

From those elements,one obtains the transitional amplitudes M^^^^, for processes J/ip —>■ pN* , N* N* by projecting 
them onto the wave functions of the proton or the Roper resonance. We ignore the trivial contribution from the 
process J/ip — » \gggg} + \gggg) to the transitional amplitude, if one assumes that Roper is the mixture of \ggg, 2s > 
state with the hybrid state \gggg >, one gets 

' ( J/V^^^ ^ pp) =< ggg, Is; W9, ls|TiCi | J/V'^''^ >, (39) 



Mf^_'^,(J/V'''' ^piV') =5<ggg,ls;gg5,2s|TiCi|J/V^^^ > +^1 - 6^ < gqg, ls;gggg\T2C2\J/i^^^'> >, (40) 



Ml'^l.iJ/i,^"^ ^ N*N*) =5< ggg,'2s;qgg,2s\TiCi\J/^'-''^ > +^1 - < ggg,2s;gggg\T2C2\J /i'^'^^ > ■ (41) 

In the decay mode J/V" — » N*N* , we only consider the process illustrated in Fig. 1(b), where the normal Roper 
state is treated as a \ggg, 2s > state, while the hybrid state is assumed to be a mixture of a 2s > state and a 
\gggg, Is > state. 

C The radial wave function for \qqq) and \qqqg) clusters 

The radial wave functions of protons or hybrid baryons are chosen as simple harmonic oscillator eigenfunctions in 
their center-of-mass (CM) system,i.e 

_> _> 1 -(fcp +kx ) 

4>is{kp,kx)= e 2„ , (42) 

(ttq:)^' 



for \qqq, Is) cluster, and 



*2 



-> -> ^1 1 ^2 ^2 -(fcp +kx ) 

,{k„kx) = V3j—^il-^ik,+kx))e ^ , (43) 



for \qqq, 2s) cluster,where a = (Sfcmq)^''^, k p, k \ are relative coordinates defined as 

= -^(^^1 + "^2 -21*3), (44) 

= (45) 

v2 

Let pi,p2,ps be momenta for three quarks and P4 for a gluon in the cluster qqqg as depicted in Fig. 3, and we define 
three relative coordinates as 

Pa = -^(pi-pi), (46) 

Pb = -^(pi +P2 - 2p3), (47) 

Pc = — ^ (pi + P2 + P3 - 3p4) . (48) 
V f2 

One has a spatial wave function for qqqg cluster in momentum space. 

_.^ 1 -iPi^+Pb^l 1 _H , 

0.(p:,P?,p2)= (^^e - __^e^., (49) 

where /3 = (SfcAfg)"'^'^^ is harmonic-oscillator parameter (Mg — )■ 




(Pl+P2)/2 



Figure 3: The definition of the relative momentum for four bodies system in CM system 
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